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Spin projection and spin current density within relativistic electronic transport 

calculations 
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A spin projection scheme is presented which allows the decomposition of the electric conductivity 
into two different spin channels within fully relativistic ab initio transport calculations that account 
for the impact of spin-orbit coupling. This is demonstrated by calculations of the spin-resolved 
conductivity of Fei-zCra; and Coi-zPtz disordered alloys on the basis of the corresponding Kubo- 
Greenwood equation implemented using the Korringa-Kohn-Rostoker coherent potential approxima- 
tion (KKR-CPA) band structure method. In addition, results for the residual resistivity of diluted 
Ni-based alloys are presented that are compared to theoretical and experimental ones that rely on 
Mott's two-current model for spin-polarized systems. The application of the scheme to deal with 
the spin-orbit induced spin Hall effect is discussed in addition. 

PACS numbers: 72.15.Eb,71.70.Ej,72.25.Ba 



During the last years research activities in spintron- 
ics increased very rapidly. The reason for the growing 
interest in this field is based on the close connection 
with fundamental scientific questions as well as its impact 
on technology [l], ■ Compared to standard electronics 
where only the charge of the electrons is used, spintronics 
uses the charge of the electrons in combination with the 
spin degree of freedom. One of the most exciting effects 
within spintronics is the spin Hall effect (SHE) 0, 
The SHE appears when an electric current flows through 
a medium with spin-orbit coupling present, leading to a 
spin current perpendicular to the charge current. This 
effect is even present in non-magnetic materials as could 
be demonstrated experimentally e.g. for Pt 
For a theoretical investigation of effects like the SHE it 
is obviously crucial to have a reliable description for the 
spin-dependent transport that accounts for the impact of 
spin-orbit coupling in a proper way. Most investigations 
in this field were based on the Pauli equation including 
spin-orbit coupling explicitly as a relativistic correction 
term and representing the spin current density essentially 
by a combination of the Pauli spin matrix a z with the 
conventional current density operator [6]. Very few in- 
vestigations have been done so far on the basis of the 
Dirac equation using an expression for the spin-current 
density, albeit introduced in an ad-hoc manner Q- In 
contrast to these approximate schemes to deal with spin- 
dependent transport the approach suggested by Vernes 
et al. [a] supplies a fully relativistic and coherent de- 
scription of electronic spin-polarization and the associ- 
ated spin-current density. This approach based on the 
four-component polarization operator T introduced by 
Bargmann and Wigner [§] leads, in particular, to a cor- 
responding set of continuity equations. 
In the present work we introduce spin projection oper- 
ators derived from the polarization operator T ■ This 
allows a decomposition of the conductivity into contri- 



butions from each spin channel within fully relativis- 
tic transport calculations. Applications on the spin- 
dependent transport of various magnetic transition metal 
alloy systems demonstrate the flexibility and reliability of 
the new approach. 

Within non-relativistic quantum mechanics the elec- 
tronic spin can be described via the well known Pauli 
matrices <7i, specifying the non-relativistic spin operator 
s = \(T. Due to the fact that the Schrodinger Hamilto- 
nian Hs commutes with s the projection of the spin, e.g. 
to the z-axis, is a constant of motion. This is no longer 
the case within a scheme that accounts for spin-orbit cou- 
pling. The most reliable approach in this context makes 
use of electronic structure calculations on the basis of 
the Dirac equation. It turns out that even in the sim- 
plest case of a free electron the Dirac Hamiltonian does 
not commute with e.g. s z . However, it is possible to de- 
fine a generalized spin operator which commutes with the 
free electron Dirac equation and shows all characteristic 
properties of a spin operator [l(| [H| ■ 
Within the fully relativistic description it is not possible 
to decompose the conductivity in a strict sense into spin- 
up and spin-down contributions in a simple way. There- 
fore, one may use approximative schemes or one can 
switch to scalar- relativistic calculations [H,[l|| to decom- 
pose the conductivity into two different spin channels. 
The disadvantage of these two approaches is that approx- 
imative schemes work only under certain circumstances 
and scalar-relativistic calculations neglect all scattering 
events that lead to a spin flip due to the fact that such 
calculations neglect spin-orbit coupling. To avoid such 
shortcomings a proper relativistic spin projection opera- 
tor is necessary. 

The starting point of our derivation of suitable relativis- 
tic spin projection operators is based on the four-vector 
polarization operator T which was derived by Bargmann 
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and Wigner [Fjj]: 
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with the kinetic momentum II = (p + — A)l4 and the 
canonical momentum p. The matrices X are the rela- 
tivistic Pauli-matrices, j3 is one of the standard Dirac 
matrices and fill: 
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The operator T can be considered as a generalized spin 
operator which commutes with the field free Dirac Hamil- 
tonian fill: 



H 



free 



cap + /3mc 
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In addition, the components 7^ are the generators of the 
little group that is a subgroup of the group of Lorentz 
transformations [l(J. In comparison to other suggested 
forms of polarization operators the operator T is gauge 
invariant [ll[ and therefore the appropriate basis for cal- 
culations which include electromagnetic fields. 
A widely used relativistic scheme to deal with magnetic 
solids within spin density functional theory was intro- 
duced by MacDonald and Vosko [T3| . The corresponding 
Dirac Hamiltonian: 



H = ca p + (3mc 2 + V cS - /3£-B cff , 



(5) 



includes an effective scalar potential V e e and an effec- 
tive magnetic field B c ff coupling only to the spin de- 
gree of freedom. For the subsequent discussion we choose 
B c ff = B{r) e z as frequently done within electronic struc- 
ture calculations. The commutator of T and H is non 
zero which shows that T is no longer a constant of mo- 
tion. 

From T corresponding spin projection operators 
can be derived by demanding: 



p++-p- = i 



or equivalently 



;(1±T) 



(0) 
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The projection of T to a unit vector along the 2-axis 
n T = (0, 0, 1, 0) leads to the following expression: 

T-n = /3£ 2 -^^. (9) 
mc 

Making use of the relation B = V x A between the 
vector potential A and the magnetic field B, A has only 
non-zero components in the xy-plane if B || e 2 (see Eq. 



([5])), i.e. A z — 0. For the spin projection operators this 
leads to: 
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Starting from the polarization operator 7~ Vernes et al. 
8] could demonstrate that a corresponding spin current 
density operator is given by a combination of T with the 
conventional relativistic electron current density opera- 
tor: 



eca b 



(11) 



where is one of the standard Dirac matrices 
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Accordingly, we get an operator for the spin-projected 
current density by combining and which leads to 



J: 



z± _ 



Using J^- to represent the observable within Kubo's 
linear response formalism one can derive expressions for 
a corresponding spin-projected conductivity tensor (the 
details will be published elsewhere). Restricting to the 
symmetric part of the tensor one arrives at: 



7 Z± 



•Tr 



J z ±^G+{E F )j v %G+{E F ) 



(12) 



Here N is the number of atomic sites, f2 the volume per 
atom, is the current density operator (p = x, y, z) and 
^G + (Ep) is the imaginary part of the retarded one par- 
ticle Green function at the Fermi energy Ep. 
Eq. (|12p is obviously the counter-part to the conventional 
Kubo-Greenwood equation 15] for the spin-integrated 
conductivity that is recovered by replacing J*^- by j^. 

For the determination of ^G + (Ep) we use multiple 
scattering theory (MST) which is the basis of the KKR 
band structure method. Within MST the real space rep- 
resentation of 3G + has the following form [l6j : 

3G+(r,r',£) = 3 £ Z% x (r n , E)t£™ A2 {E)Z^ (r m , E) , 
AiA 2 

(13) 

with r = R n + r„ and r' = R m + r m . Using a fully rela- 
tivistic implementation, the wave functions Z\ (Z^ x ) are 
the regular right (left) hand side solutions of the Dirac 
equation within cell n, t£™, is the scattering path oper- 
ator and A = (k, /1) with k and (i being the relativistic 
spin-orbit and magnetic quantum numbers (llj . 
The configurational average for a disordered alloy - indi- 
cated by the brackets (...) in Eq. (fT!2"|) - is taken by means 
of the Coherent Potential Approximation (CPA) fljj . 

The scheme outlined above has been implemented by a 
corresponding extension of the formalism worked out by 
Butler |15] to calculate the residual resistivity of disor- 
dered alloys on the basis of the Kubo-Greenwood equa- 
tion. As a first step, the underlying electronic structure 
of the investigated disordered alloy systems has been cal- 
culated self-consistently on the basis of local spin den- 
sity approximation (LSDA) using the parameterization of 
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Vosko et al. [17|. For the band structure calculations the 
fully relativistic version of the Korringa-Kohn-Rostoker 
(KKR) method [l8[ has been used in combination with 
the coherent potential approximation (CPA) alloy the- 
ory to account for chemical disorder. The CPA has been 
exploited in particular to perform the configurational av- 
erage when calculating the conductivity tensor elements 
and when dealing with the associated vertex corrections 
fl5j . These calculations have been done using a cut-off 
for the angular momentum expansion at / max = 3 to 
ensure convergence for the investigated transition metal 
systems. 

As a first application of the presented projection 
scheme the spin resolved conductivity of the alloy system 
Fei-^Cr^ has been calculated assuming the magnetiza- 
tion to be aligned along the z-axis. The presence of the 
spin-orbit coupling gives rise to the anomalous magneto 
resistance (AMR) with the conductivity tensor elements 
&xx = Cyy =/= Vzz for this situation. The reduced symme- 
try is also reflected by the spin projected conductivities 
a xx and a zz , as can be seen in Fig. [TJ Obviously, 
the conductivity is quite different for the two spin chan- 
nels. This behavior can be traced back straight forwardly 
to the electronic structure of the alloy system around 
the Fermi energy that can be represented in a most de- 
tailed way in terms of the spin-projected Bloch spectral 
function (BSF) 20]. While for the spin down subsystem 
there exists a well-defined Fermi surface with dominant 
sp-character corresponding to a sharp BSF, the spin up 
subsystem is primarily of d-character that is much more 
influenced by the chemical disorder in the system lead- 
ing to a BSF with rather washed-out features [2(J. As 
the width of the BSF can be seen as a measure for the 
inverse of the electronic lifetime the very different width 
found for the two spin subsystems explain the very dif- 
ferent spin-projected conductivities. 
Fig. [1] shows in addition results that have been obtained 
on the basis of an approximate spin-projection scheme 
that was suggested recently [Hj]. Within this scheme 
the matrices occurring in the Kubo-Greenwood equation 
for the conductivity are transformed from the standard 
relativistic representation (using the quantum numbers 
A = («, as labels) to a spin-projected one (using the 
quantum numbers L = (Z,m/,m s ) as labels). Suppress- 
ing the spin-flip term of the current density matrix el- 
ements Jlu one can easily split the conductivity into 
spin-up and spin-down contributions and an additional 
spin-flip contribution a z ^ that is related to the spin- 
off-diagonal elements of the scattering path operator r. 
For 3c?-elements with a relatively low spin-orbit coupling 
it was found that the neglect of spin-off-diagonal elements 
of Jll* is well justified and that cr zH is quite small. In 
fact the spin-projected conductivities ait and alt^ 
obtained by the approximate scheme compare very well 
with the results based on the scheme presented here (see 
Fig. EE]). 

In order to demonstrate the limitations of the approx- 
imate scheme from Ref. [lj| the isotropic spin resolved 
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FIG. 1: (Color online) Spin resolved conductivity tensor el- 
ements cr xx and a zz of Fei-^Crz calculated for the 
magnetization pointing along the z-axis (full symbols). In 
addition, results are shown that have been obtained using an 
approximate scheme (open symbols) [T^ ]. 



conductivity for Coi-^Pta, is shown in Fig. [5J It turns 
out, that the approximative scheme fails especially for 
high Pt concentrations. This can be attributed to an in- 
creased spin-orbit interaction for which the assumptions 
on which this scheme is based are no longer fulfilled. 

As another application of the presented scheme results 
for diluted Ni-based alloys with XNi = 0.99 are shown 
in Fig. [3] in terms of the isotropic residual resistivities 



P 



*+(-) 



[(2 0". 



*+(- 

XX 



,*+(-) 



)/3]" 



As one notes, the 



resistivity for the two spin channels show a rather differ- 
ent variation with the atomic number of the impurities. 
This can be traced back again to the spin-projected elec- 
tronic structure of Ni at the Fermi level and the position 
of the impurity d-states 12 111 . In Fig. [3] the results of cal- 
culations by Mertig et al. [21| have been added, that were 
done in a scalar-relativistic way - i.e. ignoring spin-orbit 
coupling - on the basis of the Boltzmann-formalism and 
by making use of the two-current model. In spite of the 
various differences between this approach and the pre- 



4 



200 
160 



80 



40 



1 1 

- ^ o-o a total 
\ A - A a (approx.) 
\ v^a (approx.) 
\ A *G Z+ (proj.) 
\ a z " (proj.) J 


1 




:-r— 

1 



0.2 



0.4 0.6 

x 



0.8 



FIG. 2: (Color online) Isotropic spin resolved conductivity 



,*+(-) 



(2(j^' _ ^ + crj^ 1 - j )/3 of Coi-zPt^ for the magne- 
tization pointing along the z-axis (full symbols). In addition, 
results are shown that have been obtained using an approxi- 
mate scheme (open symbols) [l9| . 



scntcd scheme, the resulting spin-projected resistivities 
agree fairly well. This also holds concerning correspond- 
ing experimental data that have been deduced from mea- 
surements relying on the two-current model. 

In summary, a scheme for a spin projection within 
transport calculations on the basis of the Kubo formal- 
ism has been presented. The applications presented were 
restricted to the diagonal elements of the corresponding 
conductivity tensor described by a Kubo-Grccnwood-likc 
equation. Results obtained for the disordered alloy sys- 
tems Fei-^Cra;, Coi^Ptj. and diluted Ni-based alloys 
were compared to results based on an alternative but 
approximate projection scheme and theoretical as well 
as experimental data based on the two current model. 
The good agreement found for the investigated systems 
ensures the consistency and reliability of the presented 
scheme. Accordingly, this is expected to hold also when 
dealing with spin-projected off-diagonal conductivities as 

z+( — ) 

e.g. a xy on the basis of Kubo-Stfeda-like equations. 
This will give access in particular to the spin-projected 
Hall conductivity in magnetic materials as well as to the 
spin Hall conductivity in non-magnetic materials. Work 
along this line is in progress. 
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FIG. 3: (Color online) Isotropic spin resolved resistivity of 
Ni with 3d transition metal impurities (1%) obtained by the 
present scheme (full squares) compared to theoretical data 
from Mertig et al. I2E (blue squares/dotted line), experimen- 
tal data from Fert |22| (full red circles/dashed line) and other 
experimental data (see Ref. |23| . open red circles). The top 
and the bottom panel show the data for spin-up and for spin- 
down, respectively. 
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